Spin-pumping and Enhanced Gilbert Damping in Thin Magnetic Insulator Films 
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Precessing magnetization in a thin film magnetic insulator pumps spins into adjacent metals; 
however, this phenomenon is not quantitatively understood. We present a theory for the dependence 
of spin-pumping on the transverse mode number and in-plane wave vector. For long-wavelength spin 
waves, the enhanced Gilbert damping for the transverse mode volume waves is twice that of the 
macrospin mode, and for surface modes, the enhancement can be ten or more times stronger. Spin- 
pumping is negligible for short-wavelength exchange spin waves. We corroborate our analytical 
theory with numerical calculations in agreement with recent experimental results. 
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Metallic spintronics have been tremendously success- 
ful in creating devices that both fulfill significant market 
needs and challenge our understanding of spin transport 
in materials. Topics that are currently of great interest 
are spin transfer and spin-pumping [lH3[, spin Hall ef- 
fects [4[ , and combinations thereof for use in non- volatile 
memory, oscillator circuits, and spin wave logic devices. 
A recent experimental demonstration that spin transfer 
and spin-pumping can be as effective in magnetic insula- 
tors as in metallic ferromagnetic systems was surprising 
and has initiated a new field of inquiry [f| . 

In magnetic insulators, no moving charges are present, 
and in some cases, the dissipative losses associated with 
the magnetization dynamics are exceptionally low. Nev- 
ertheless, when a magnetic insulator is placed in con- 
tact with a normal metal, magnetization dynamics in- 
duce spin-pumping, which in turn causes angular momen- 
tum to be dumped to the metal's itinerant electron sys- 
tem. Due to this non-local interaction, the magnetization 
losses become enhanced. Careful experimental investiga- 
tions of spin-pumping and the associated enhanced mag- 
netization dissipation were recently performed, demon- 
strating that the dynamic coupling between the magne- 
tization dynamics in magnetic insulators and spin cur- 
rents in adjacent normal metals is strong. Importantly, in 
magnetic insulators, an exceptionally low intrinsic damp- 
ing combined with good material control has enabled the 
study of spin-pumping for a much larger range of wave 
vectors than has previously been obtained in metallic fer- 
romagnets foUl^. 

In thin film ferromagnets, the magnetization dynamics 
are strongly affected by the long-range dipolar interac- 
tion, which has both static and spatiotemporal contribu- 
tions. This yields different types of spin waves. When 
the in-plane wavelength is comparable to the film thick- 
ness or greater, the long-range dipolar interaction causes 
the separation of the spin-wave modes into three classes 
depending on the relative orientation of the applied ex- 
ternal field, in relation to the film normal, and the spin- 
wave propagation direction [lil - Eoj . These spin waves 
are classified according to their dispersion and transverse 
magnetization distribution as forward volume magneto- 



static spin waves (FVMSWs) when the external field is 
out-of-plane, backward volume magnetostatic spin waves 
(BVMSWs) when the external field is in-plane and along 
the direction of propagation, and magnetostatic surface 
waves (MSSWs) when the external field is in-plane but 
perpendicular to the direction of propagation. In volume 
waves, the magnetic excitation is spatially distributed 
across the entire film, while surface modes are localized 
near one of the surfaces. "Backward" waves have a fre- 
quency dispersion with a negative group velocity for some 
wavelengths. While these spin waves have been studied 
in great detail over the last decades, the effect of an adja- 
cent normal metal on these waves has only recently been 
investigated. 

Experimentally, it has been observed that spin- 
pumping differs for FVMSWs, BVMSWs and MSSWs 
and that it depends on the spin- wave wavelength [f| [|| 0, 
12 "14 1- Recent experiments [1] have demonstrated that 
the magnetization dissipation is larger for surface spin 
waves in which the excitation amplitude is localized at 
the magnetic insulator-normal metal interface. To uti- 
lize spin-pumping from thin film magnetic insulators into 
adjacent normal metals, a coherent theoretical picture of 
these experimental findings must be developed and un- 
derstood, which is the aim of our work. 

In this Letter, we present a theory for energy dissipa- 
tion from spin- wave excitations in a ferromagnetic insula- 
tor (FI) thin film via spin-pumping when the ferromag- 
netic insulator layer is in contact with a normal metal 
(NM). To this end, consider a thin film magnetic insu- 
lator of thickness L on an insulating substrate with a 
normal metal capping (see Figure Q]). We consider a nor- 
mal metal such as Pt at equilibrium, where there is rapid 
spin relaxation and no back-flow of spin currents to the 
magnetic insulator. The normal metal is then a perfect 
spin sink and remains in equilibrium even though spins 
are pumped into it. 

The magnetization dynamics are described by the 
Landau-Lifshitz-Gilbert (LLG) equation [2l[ with a 
torque originating from the FI/NM interfacial spin- 
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where a is the Gilbert damping coefficient, Ms is the 
saturation magnetization, 7 is the gyromagnetic ratio, 
H e ff is the effective field including the external field, ex- 
change energy, surface anisotropy energy, and static and 
dynamic demagnetization fields. 

Spin-pumping through interfaces between magnetic in- 
sulators and normal metals gives rise to a spin-pumping 
induced torque that is described as @ 
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where g± is the transverse spin ( "mixing" ) conductance 
per unit area at the FI/NM interface. We disregard the 
imaginary part of the mixing conductance because this 
par t has been found to be small at FI/NM interfaces 
12j . In addition, the imaginary part is qualitatively less 
important and only renormalizes the gyromagnetic ratio. 

Assuming only uniform magnetic excitations, 
"macrospin" excitations, the effect of spin-pumping 
on the magnetization dissipation is well known [2 
Spin-pumping leads to enhanced Gilbert damping, 
a — > a + Aa macro , which is proportional to the FI/NM 
cross section because more spin current is then pumped 
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out, but inversely proportional to the volume of the 
fcrromagnet that controls the total magnetic moment: 
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Thus, the enhanced Gilbert damping due to spin- 
pumping is inversely proportional to the film thickness 
L and is important for thin film ferromagnets. However, 
a "macrospin" excitation, or the FMR mode, is only one 
out of many types of magnetic excitations in thin films. 
The effect of spin-pumping on the other modes is not 
known, and we provide the first analytical results for this 
important question, which is further supported and com- 
plemented by numerical calculations. 

We consider weak magnetic excitations around a ho- 
mogenous magnetic ground state pointing along the di- 
rection of the internal field = if;Z, which is the com- 
bination of the external applied field and the static de- 
magnetizing field [3] . We may then expand M = Msz + 
m Q,2y(£) eiM " QC) > where m Q , xy ■ z = 0, \m QlXy \ < Ms, 
and Q is the in-plane wave number in the ^-direction. 

Following the linearization approach of the LLG equa- 
tion (JTJ as in Ref. [l9j], we arrive at a two-dimensional 
intcgro-differcntial equation of the dynamic magnetiza- 
tion (in the xy-plane) in the film's transverse coordinate 
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where lo is the spin-wave eigenfrequency, A is the ex- 
change stiffness, ujh = 7 Hi, ujm = AirjMs, and Q xy is 
the dipole-dipole field interaction tensor, which fulfills 
the boundary conditions resulting from Maxwell's equa- 
tions (see [23j). 

The cigensystem must be supplemented by boundary 
conditions that account for spin-pumping and surface 
anisotropy. These boundary conditions are obtained by 
integrating Eq. ([T]) over the interface 24| and expanding 
to the lowest order in the dynamic magnetization. When 
an out-of-plane easy axis surface anisotropy is present, 
the boundary conditions are 
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where K s is the surface anisotropy energy with units 
erg ■ cm -2 and \ = ^ s a parameter relating the ex- 



change stiffness and the spin mixing conductance ([x] = 
s). The boundary condition at the magnetic insulator- 
substrate interface at £ = —L/2 is similar to Eq. ([5]), but 
simpler because x and K s — » at that interface. 

A mathematical challenge induced by spin-pumping 
arises because the second term in the linearized boundary 
condition ([5]) is proportional to the eigenvalue u such that 
the eigenfunctions cannot simply be expanded in the set 
of eigenfunctions obtained when there is no spin- pumping 
or dipolar interaction. Instead, we follow an alternative 
analytical route for small and large wave vectors. Fur- 
thermore, we numerically determine the eigenmodes with 
a custom-tailored technique, where we discretize the dif- 
ferential equation Q, include the spin-pumping bound- 
ary conditions (O, and transform the resulting equations 
into an eigenvalue problem in w [2o| . 

Let us now outline how we obtain analytical results for 
small QL <C 1 and large QL 1 wave vectors. First, 
we consider the case of vanishing surface anisotropy and 
compute the renormalization of the Gilbert damping for 
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(see Eq. ©), Aa = Aa„ 
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FIG. 1. a) A thin film magnetic insulator of thickness L 
in its coordinate system; £ is the normal axis, the infinite 
77^-plane is coplanar with the interfaces, and the spin waves 
propagate along the ("-axis. The internal field and saturation 
magnetization are along the z-axis. The y-axis is always kept 
in-plane, and the a;-axis is selected such that the x-, y- and 
z-axes form a right-handed coordinate system, b) A cross- 
section showing the material stack. 



the resulting modes. Next, we demonstrate that the sur- 
face anisotropy creates a surface wave with a comparably 
large enhancement of the Gilbert component. 

When QL <C 1, the convolution integral on the right- 
hand side of Eq. (TJJ only contains the homogeneous de- 
magnetization field. The magnetization is then a trans- 



verse standing wave m.Q^ xy (e 



, where k is 



a transverse wave number, <fi is a phase determined by 
the BC at the lower interface, and the two-dimensional 
coefficient vector m.Q^ xy allows for elliptical polarization 
in the xy-plane. 

By employing exchange-only boundary conditions [24j 
at the lower interface and using Eq. ([5]) with K s = on 
the upper interface, the transverse wave number k is de- 
termined by kLtankL — iojx- Together with the bulk 
dispersion relation uo = uj(k), calculated from Eq. Q, 
this expression allows us to calculate the magnetic excita- 
tion dispersion relation parameterized by the film thick- 
ness, the Gilbert damping a, and the transverse conduc- 
tance g±. 

When spin-pumping is weak, cox is small, and the so- 
lutions of the transcendental equation can be expanded 
around the solutions obtained when there is no spin- 
pumping, kL = nir, where n is an integer. When 
n 0, we expand to first order in kL and obtain 
kL « nir + icox/(nn). When n = 0, we must perform 
a second-order expansion in terms of kL around 0, which 
results in (kL) 2 w i^X- Using these relations in turn 
to eliminate k from the bulk dispersion relation while 
maintaining our linear approximation in small terms and 
solving for cj, we obtain complex eigenvalues, where the 
imaginary part is proportional to a renormalized Gilbert 
damping parameter, a* = a + Aa. When n = 0, our 
results agree with the spin-pumping-enhanced Gilbert 
damping of the macrospin (FMR) mode derived in [2| 



These new results indicate that all higher transverse vol- 
ume modes have an enhanced magnetization dissipation 
that is twice that of the macrospin mode. Thus, coun- 
tcrintuitivcly, with the exception of the macrospin mode, 
increasingly higher-order standing-wave transverse spin- 
wave modes have precisely the same enhanced Gilbert 
damping. 

Next, let us discuss spin-pumping for surface waves 
induced by the presence of surface anisotropy. When 
K s =^= 0, the lowest volume excitation mode develops into 
a spatially localized surface wave. Expanding the ex- 
pression for the localized wave to the highest order in 
LK S /A, we determine after some algebra that the result- 
ing enhancement of the Gilbert damping is 



COM 



K: 



4nM 2 A 



jhK s h co H 
Aan=0 ~ 4^KA^ g± ^ 

(7) 

Comparing Eqs. (O and ©, we see that for large sur- 
face anisotropy LK S /A 3> 1, the spin-pumping- induced 
enhanced Gilbert damping is independent of L. This re- 
sult occurs because a large surface anisotropy induces 
a surface wave with a decay length A/K s , which re- 
places the actual physical thickness L as the effective 
thickness of the magnetic excitations, i.e., for surface 
waves L — > A/K s in the expression for the enhanced 
Gilbert damping of Eq. (|3|). This replacement implies 
that the enhanced Gilbert damping is much larger for 
surface waves because the effective magnetic volume de- 
creases. For typical values of A and K s , we obtain an 
effective length A/K s ~ 10 nm. Compared with the film 
thicknesses used in recent experiments, this value corre- 
sponds to a tenfold or greater increase in the enhance- 
ment of the Gilbert damping. In contrast, for the volume 
modes (n / 0), we note from Eq. ((5|) that the dynamic 
magnetization will decrease at the FI/NM interface due 
to the surface anisotropy; hence, Aa decreases compared 
with the results of Eq. ([6]) . 

Finally, we can also demonstrate that for large wave 
vectors QL 3> 1, the excitation energy mostly arises from 
the in-plane (longitudinal) magnetization texture gradi- 
ent. Consequently, spin-pumping, which pumps energy 
out of the magnetic system due to the transverse gradi- 
ent of the magnetization texture, is much less effective 
and decays as 1/(QL) 2 with respect to Eq. J3]). 

To complement our analytical study, we numerically 
computed the eigenfrequencies ui n (Q). The energy is de- 
termined by the real part of uj n (Q), while 3mu n (Q) de- 
termines the dissipation rate and hence the sp in-pump ing 
contribution. Recent experiments 0, [HI [HI, Il4l | on 
controlling and optimizing the ferrimagnetic insulator 
yttrium-iron-garnet (YIG) have estimated that the mix- 
ing conductances of both YIG — Au and YIG — Pt bilayers 
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FIG. 2. Aa versus wave vector for the MSSW geometry (6 = 
(j> = 7r/2) for the four lowest eigenvalues. Inset: Magnitudes of 
eigenvectors (in arbitrary units) across the film at QL — 1.5. 



are in the range of g±h/e 2 ~ 0.02-3.43 • 10 15 cm" 2 . We 
use g±h/e 2 = 1.2 • 10 14 cm~ 2 from Ref. @ in this work. 
All of our results can be linearly re-scaled with other val- 
ues of the mixing conductance. In the following section, 
we also use A = 2.9 • 10 -8 erg/cm, K s = 0.05 erg/cm 2 , 
L = 100 nm, 4ttA/ s = 1750 G, and a = 3 ■ 10" 4 . 

To distinguish the spin-pumping contribution Aa from 
the magnetization dissipation due to intrinsic Gilbert 
damping a, we first compute the eigenvalues, lj^, with 
intrinsic Gilbert damping, and no spin-pumping, 

g± = 0. Second, we compute the eigenvalues w sp with 
dissipation arising from spin-pumping only, a = and 
g± =/= 0. Because 3m oc a, we define a measure of 
the spin-pumping-induccd effective Gilbert damping as 
Aa = a 3mu! sp / 3mwd- 

We first consider the case of no surface anisotropy. Fig- 
ure[5]shows the spin-pumping-enhanced Gilbert damping 
Aa as a function of the product of the in-plane wave vec- 
tor and the film thickness QL in the MSSW geometry. 
In the long- wavelength limit, QL <C 1, the numerical re- 
sult agrees precisely with our analytical results of Eq. (|6|) . 
The enhanced Gilbert damping of all higher transverse 
modes is exactly twice that of the macrospin mode. In 
the dipole-exchange regime, for intermediate values of 
QL, the dipolar interaction causes a small asymmetry in 
the eigenvectors for positive and negative eigenfrcqucn- 
cies because modes traveling in opposite directions have 
different magnitudes of precession near the FI/NM in- 
terface [26|, and spin- pumping from these modes there- 
fore differ. This phenomenon also explains why the en- 
hanced damping, Aa, splits into different branches in 
this regime, as shown in Fig. [2] For exchange spin waves, 
QL 3> 1, the exchange interaction dominates the dipo- 
lar interaction and removes mode asymmetries. We also 
see that Aa — > for large QL, in accordance with our 
analytical theory. 

Figure [3] shows Aa for the BVMSW geometry. The 
eight first modes are presented; however, as no substan- 
tial asymmetry exists between eigenmodes traveling in 



FIG. 3. Aa versus wave vector for the BVMSW geometry 
(9 = 7r/2 and <j> — 0). Left inset: Magnitude of eigenvectors 
(in arbitrary units) across the film when QL — 1.5. Right 
inset: The real part of the dispersion relation for the same 
modes. 
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FIG. 4. Aa versus wave vector for the MSSW geometry 
(8 = 4> = f/2) with surface anisotropy added at the inter- 
face. Inset: Magnitudes of eigenvectors (in arbitrary units) 
across the film. 



different directions, the modes have the same pairwise 
renormalization of a. This symmetry occurs because the 
direction of the internal field coincides with the direction 
of propagation. As in the previous case, the dipolar in- 
teraction causes a slight shift in the eigenvectors in the 
intermediate QL regime, thereby altering Aa from that 
of Eq. @. 

Figure S] shows Aa for the MSSW geometry but with 
surface anisotropy at the FI/NM interface. As expected 
from our analytical results, surface anisotropy induces 
two localized surface modes with a ten-fold larger en- 
hancement of Aa compared with the volume modes. The 
horizontal dashed line in Figure 0] indicates the analyti- 
cal result for the enhanced Gilbert damping of the ti/0 
modes when K s = 0. For the volume modes, it is clear 
that the eigenvectors have a lower magnitude closer to the 
FI/NM interface and that Aa is lower compared with the 
case of K s = 0, which is consistent with our analytical 
analysis. 

Our results also agree with recent experiments. 
Sandweg et al. found that spin-pumping is signifi- 
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cantly higher for surface spin waves compared with vol- 
ume spin- wave modes. In addition, in Ref. @, exchange 
waves were observed to be less efficient at pumping spins 
than dipolar spin waves, which is consistent with our re- 
sults. Furthermore, our results are consistent with the 
theoretical finding that spin-transfer torques preferen- 
tially excite surface spin waves with a critical current 
inversely proportional to the penetration depth (27l |. 

In conclusion, we have analyzed how spin-pumping 
causes a wave-vector-dcpcndcnt enhancement of the 
Gilbert damping in thin magnetic insulators in con- 
tact with normal metals. In the long-wavelength limit, 
our analytical results demonstrate that the enhancement 
of the Gilbert damping for all higher-order volumetric 
modes is twice as large as that of a macrospin excita- 
tion. Importantly, surface anisotropy-pinned modes have 
a Gilbert renormalization that is significantly and lin- 
early enhanced by the ratio LK 8 /A. 
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